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In our previous paper [Phys. Rev. A 84 042303 (2011)], we proposed an efficient scheme to 
purify dynamically a bipartite entangled state using short chains of atoms coupled to high-finesse 
optical cavities. In contrast to conventional entanglement purification protocols, we avoid controlled- 
NOT gates and thus reduce complicated pulse sequences and superfluous qubit operations. In this 
paper, we significantly improve the output fidelity of remotely entangled atoms by introducing one 
additional entanglement protocol in each of the repeater nodes and by optimizing the laser beams 
required to control the entire scheme. Our improved distillation scheme yields an almost unit 
output fidelity that, together with the entanglement distribution and swapping, opens an attractive 
route towards an efficient and experimentally feasible quantum repeater for long-distance quantum 
communication. 
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I. INTRODUCTION 

In classical data transmission, repeaters are used to 
amplify the data signals (bits) when they become weaker 
during their propagation. In contrast to classical infor- 
mation, the above mechanism is impossible to realize 
when the transmitted data signals carry bits of quan- 
tum information (qubits). In an optical- fiber system, for 
instance, a qubit is typically encoded by a single pho- 
ton which cannot be amplified or cloned without destroy- 
ing quantum coherence associated with this qubit [TJ |2] . 
Therefore, the photon has to propagate along the entire 
length of the fiber, which causes an exponentially de- 
creasing probability to receive this photon at the end of 
the channel. 

To avoid this exponential decay of a photon wave- 
packet and preserve its quantum coherence, the quan- 
tum repeater was proposed [3 . This repeater can be 
divided in three building-blocks which have to be ap- 
plied sequentially. First, a large set of entangled photon 
pairs distributed over sufficiently short fiber segments are 
generated. The two subsequent steps, (i) entanglement 
purification [H [5] and (ii) entanglement swapping [6 , are 
employed to extend the short-distance entangled photon 
pairs over the entire length of the channel. Using the 
entanglement purification, high-fidelity entangled pairs 
are distilled from a larger set of low-fidelity entangled 
pairs by means of local operations performed in each of 
the repeater nodes and classical communication between 
these nodes. The entanglement swapping, finally, com- 
bines two entangled pairs of neighboring segments into 
one entangled pair, gradually increasing the distance of 
shared entanglement. 
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Because of the fragile nature of quantum correlations 
and inevitable photon loss in the transmission channel, in 
practice, it poses a serious challenge to outperform the di- 
rect transmission of photons along the fiber. Up to now, 
only particular building blocks of a quantum repeater 
have been experimentally demonstrated, i.e., bipartite 
entanglement purification ^ ^ , entanglement swapping 
[9l|l0], and entanglement distribution between two neigh- 
boring nodes [lT[ |T2] . Motivated both by an impressive 
experimental progress and theoretical advances, more- 
over, various revised and improved implementations of 
repeaters and their building-blocks have been recently 
proposed p^HTTj . 

Practical schemes for implementing a quantum re- 
peater are not straightforward. The two mentioned 
protocols, entanglement purification and entanglement 
swapping, in general, require feasible and reliable quan- 
tum logic, such as single- and two-qubit gates. Because 
of the high complexity and demand of physical resources, 
entanglement purification is the most delicate and cum- 
bersome part of a quantum repeater. The conventional 
purification protocols [5l[T8], moreover, involve multiple 
controlled-NOT (CNOT) gates which pose a serious chal- 
lenge for most physical realizations of qubits, involving 
complicated pulse sequences and superfluous qubit oper- 
ations [HiiMia]. 

In our previous paper [24], we suggested a more prac- 
tical scheme to purify a bipartite entangled state by ex- 
ploiting the natural evolution of spin chains instead of 
CNOT gates. The realization of this dynamical scheme 
was proposed in the framework of cavity QED using short 
chains of atoms and optical cavities. In the present pa- 
per, we propose a modified purification scheme, in which 
we significantly improve the output fidelity of remotely 
entangled atoms. By introducing one additional entan- 
glement protocol in each repeater node and by optimizing 
the laser beams required to control the entire scheme, we 
reach an almost unit output fidelity after the same num- 
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FIG. 1. (Color online) (a) Sequence of steps in the original purification scheme, (b) Quantum circuit corresponding to the 
interaction indicated above by grey ellipses, (c) Structure of a three-level atom in the A- type configuration, (d) Experimental 
setup that realizes the purification scheme (a) and is incorporated into a quantum repeater segment with two neighboring 
nodes. See text for description. 



ber of purification rounds as before. This dramatic im- 
provement, therefore, allows for multiple entanglement 
swapping operations and opens a route towards an effi- 
cient and experimentally feasible quantum repeater for 
long-distance quantum communication. 

The paper is organized as follows. In the next section, 
we describe in detail the original purification scheme pre- 
sented in our previous paper. In Sec. Ill, we present our 
modified high-fidelity purification scheme. We analyze 
the atomic evolution mediated by the cavity and laser 
field, and we determine the main properties which are 
relevant for our scheme in Sec. III. A. In Sees. III.B and 
IV, we discuss a few relevant issues related to the im- 
plementation of our purification scheme, while a short 
summary and outlook are given in Sec. V. 



II. DYNAMICAL ENTANGLEMENT 
PURIFICATION 

In our previously proposed purification scheme, two 
repeater nodes A and B share one permanent qubit pair 



Pab and a finite set of temporary (low-fidelity) entangled 
pairs grouped into elementary blocks of two qubit pairs 
as displayed in Fig.jlja). Each temporary entangled pair 
is given by the rank- two mixed state 



Pf 



(1) 



where ^ 



\^ab){^ab\ -^^^^ states in the 

qubit- storage basis {|0), |1)}, and where the fidelity 



F(p/)=Tr ^X^pf =/>0.5 



(2) 



is above the threshold value of 1/2. The qubit-storage 
states |0) and |1) correspond to the two long-living states 



^ In our previous paper, we considered the Werner state pf = 
+ ^^^(<l>~ + ^+ + ^~) to describe low-fidelity entangled 
pairs. In this paper, instead, we consider the state ([T]) that can 
be efficiently generated between two remote nodes of a repeater 
using an optimal, ultimate entanglement distribution and detec- 
tion protocol ^|. In the last section, we discuss this protocol 
and provide evidence supporting our choice. 
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FIG. 2. (Color online) (a) Fidelities /) (solid curve) and 
F3{f,F2) (dashed curve) given by Eqs. (|8| and (11), respec- 
tively. The dotted curve displays the fidelity give n by Eq. (34) 
in Ref. [24]. (b) Plot of F{f, n) given by Eq. ([lo]) as a function 
of input fidelity / and number n of purification rounds. 



qubits. The above Hamiltonian with periodic bound- 
aries is produced deterministically in our scheme by cou- 
pling of three (three-level) atoms to the same mode of 
a high-finesse resonator [see Fig. [iJcY]. In our previous 
paper, we identified the Hamiltonian (p| with the Jaynes- 
Cummings Hamiltonian in the large detuning limit, i.e., 
A ^ where g is the atom-cavity coupling strength, 
A is the atom-cavity detuning, and Ji = ^^/A is the 
coupling between the atoms subject to the same cavity 
mode. 

The evolution governed by the Hamiltonian ([3| over 
the time period Q is referred to below as the purifi- 
cation gate and is indicated by the ellipses in Fig. [TJa) 
and by rectangles in Fig. [ijb). Once the purification 
gate is performed, the qubit pairs Iab and 2ab are 
pairwise projected in the computational (qubit-storage) 
basis {|0), |1)} and the outcome of these projections is 
exchanged between the two nodes by means of classi- 
cal communication [see Fig. [ijb) and the third box of 
Fig. [TJa)]. Entanglement purification is successful if the 
outcome of projections reads 

{0ia,12A,0ib,12b} or {1ia, O2A, liB, 02b} • (5) 

In this case, the (unprojected) permanent qubit pair Pab 
is described by the density operator 

PF = Fif, n $+ ^ + (1 - Fif, /')) , (6) 

where 



of a three-level atom in the A-configuration as displayed 
in Fig. [ijc). In order to protect this qubit against the 
decoherence caused by the fast-decaying excited state |e), 
the states |0) and |1) are chosen as the stable ground 
and long-living metastable states or as the two hyperfine 
levels of the ground state. 

The permanent pair Pab , characterized by the density 
operator pj/, is supplemented by two temporary pairs 
Iab and 2ab^ characterized by the density operators pj- 
and pj, respectively, as seen Fig. [ija). Each of the re- 
peater nodes A and 5, therefore, contains one triplet of 
qubits Pa, 1^,2^ and Pb^Ib^'^b^ respectively. Each of 
these triplets evolves due to the isotropic Heisenberg XY 
Hamiltonian [26l 



i=l 



over the time period (n = 0, 1, 2, . . .) 



(3) 



(4) 



where af and erf are the respective Pauli operators in 
the cavity- active basis {|0), |e)}, such that af = and 
cf\ = cr^, and where Ji is the coupling between the 



Fif,n = 



f -16(f -2)/ + 32(3f -1)/^ 
81 + 32/2 _ 80/' + 16(10/' - 7)/ ' 



(7) 



such that F(pi?) > F(pj/). The entanglement purifica- 
tion is unsuccessful, if the mentioned outcome of pro- 
jections Iab and 2ab disagrees with ([5|. In this case, 
the permanent pair Pab should be reinitialized and the 
entire sequence from Fi g, p^ a) restarted. 

The density operator ([6|) ensures that the (permanent) 
qubit pair Pab preserves its rank-two form after each suc- 
cessful purification round. Unlike the conventional purifi- 
cation protocol, therefore, the purified state ^ is com- 
pletely characterized by the fidelity F (pp) = F{f^f'). 
The expression ([7|), furthermore, describes quantitatively 
how the input fidelity of the permanent qubit pair is 
modified due to one single (and successful) purification 
round. In Fig. [2|a), we compare the fidelity 



F{fJ) 



/(II -16/ + 32/2) 
27 - 64/ + 64/2 



(8) 



(solid curve) with the respective fidelity given by Eq. (34) 
in Ref. [24 (dotted curve) that was obtained within the 
same scheme, however, by considering the Werner state 
instead of the rank- two mixed state ([T]) in this paper. As 
expected due to the vanishing contribution of in ([l]), 
the growth of fidelity in the case of a rank-two mixed 
state is larger as for the Werner state. 
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Assuming that each purification round is successful, 
the sequence from Fig. [TJa) leads to the gradual growth 
of entanglement fidelity (of stationary atoms) with regard 
to the respective fidelity obtained in the previous round 

/' < /') < F2(/, Fi) < . . . < F„(/, F„_i). (9) 

In order to understand how much the output fidelity in- 
creases with each purification round, we analyze quanti- 
tatively the following sequence 

/< /)<... < i^n(/, Fn-i) = / + n) . (10) 



In Fig. jij^b), we show the plot of function F{f^n) 
that describes the difference between the final fidelity 
Fn{f^ Fn-i) obtained after n (successful) rounds and the 
initial fidelity / (n = 0). It is clearly seen that during 
the first three rounds, this function exhibits a notably 
fast growth that saturates and, with increasing n, yields 
a negligible growth with regard to the fixed point fidelity 



I 

/ (70859 - 377904/ + 950112/^ - 1368064/^ + 1278976/"^ - 671744/^ + 294912/^) 
177147 - 1051072/ + 2792896/^ - 4204544/^ + 3904512/4 - 2162688/^ + 720896/^ ' 



Regardless of the number of (successful) purification 
rounds, therefore, the final fidelity is bounded by a fixed 
point value that determines the optimal number of purifi- 
cation rounds required to reach the best performance in 
a resource- and time-efficient way. Such behavior is the 
common feature of the (so-called) entanglement pump- 
ing purification scheme introduced by W. Diir and co- 
authors in Ref. [18 . We refer to this property as the 
saturation of entanglement purification. Corresponding 
to n = 3 successful purification rounds, for which the 
final fidelity reaches its saturation level, we display in 
Fig. [2|a) the (fixed point) fidelity F3(/, F2) by a dashed 
curve. By comparing this curve to the solid curve, we 
conclude that the growth of fidelity due to three succes- 
sive rounds is notably larger if compared to the case of a 
single purification round. 

We remark that the described purification scheme is 
based on the effect of entanglement transfer between the 
networks of evolving spin chains that was introduced and 
investigated in Ref. [27 . In the same reference, it was 
suggested that this effect plays the key role in the en- 
tanglement concentration once a part of the spins from 
two such networks are locally measured. In our scheme, 
the role of (spin-chain) networks is played by the atomic 
triplets located in two repeater nodes, while the cavity- 
mediated interaction governed by the Hamiltonian (|3| 
reproduces the spin-chain dynamics. From a more fun- 
damental point of view, the mentioned effect of entangle- 
ment transfer originates the constructive and destructive 
interference of the quantized spin waves (magnons) in an 
evolving spin chain (see [22 and references therein). 

The main physical resources of the proposed purifica- 
tion scheme are: (i) short chains of atoms, (ii) two high- 
finesse optical cavities, and (iii) detectors for projective 
measurements of atomic states. In Fig. [ijd) we show 
the experimental setup of a quantum repeater segment 
that includes two neighboring nodes (A and B). In this 
setup, each repeater node consists of one optical cavity 
Ci {C2) acting along the ?/-axis, a laser beam Li (I/3), a 



chain of atoms transported by means of an optical lattice 
along the same axis, one stationary atom trapped inside 
the cavity with the help of a vertical lattice, laser beam 
L2 {L4) acting along the ?/-axis, a magneto-optical trap 
(MOT), and a CCD camera connected to the neighboring 
node through a classical communication channel. 

We associated the permanent qubits with the station- 
ary atoms trapped inside cavities Ci and C2, and the 
temporary qubits with (the chains of) atoms inserted 
into the horizontal lattices and transported along the z- 
axis. According to the experimental scheme in Fig. [TJd), 
this identification implies that atoms pass sequentially 
through the cavity, such that only two atoms from the 
chain couple simultaneously to the same cavity mode. 
These two atoms together with the stationary (trapped) 
atom form an atomic triplet in each repeater node as 
assumed by our purification scheme. 

Right before an atom from node A enters the cavity, it 
becomes entangled with the respective atom from node B 
as depicted in Fig. TJd) by wavy lines. This entanglement 
is generated non-locally by means of an entanglement 
distribution block (indicated by a rectangle), such that 
each produced entangled pair is described by Eq. ([T]) in 
the qubit-storage basis {|0), |1)}. During the transition of 
an atomic pair through the cavity, the triplet of atoms has 
to undergo the cavity-mediated evolution governed by the 
Hamiltonian pi) in each of the repeater nodes over the 
time period Since the Hamiltonian acts solely on the 
cavity- active states {|0), |e)}, the atomic population has 
to be mapped from the qubit-storage basis to the cavity- 
active basis in order to make possible the interaction of 
atoms with the cavity mode and, moreover, to protect 
the qubits against the decoherence caused by the fast- 
decaying excited state |e). This mapping is realized using 
short resonant light pulses produced by the laser beam Li 
(L3). Each pulse transfers the electronic population from 
the qubit-storage states to the cavity-active states (or 
backwards), such that the atoms couple to (or decouple 
from) the cavity field in a controlled fashion. 
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FIG. 3. (Color online) (Color online) (a) Sequence of steps in the modified purification scheme, (b), (c) Structure of a three-level 
atom in the A-type configuration subjected to the cavity and laser fields. See text for description, (d) Experimental setup that 
realizes the purification scheme (a) and is incorporated into a quantum repeater segment with two neighboring nodes. 



According to the sequence in Fig. TJa), furthermore, 
the purification sequence is completed once the states of 
an (conveyed) atomic pair are projectively measured and 
the outcome of projections is pairwise exchanged between 
the repeater nodes in order to decide if the purification 
was successful or not. In our experimental scheme, the 
latter projections are performed by means of the laser 
beam L2 (I/4) and a CCD camera in each of the repeater 
nodes as displayed in Fig. [ijd). While the laser beam 
L2 {L4) removes atoms in a given (storage-basis) state 
from the chain without affecting atoms in the other state 
(so-called push-out technique [2F), the CCD camera is 
used to detect the presence of remaining atoms via flu- 
orescence imaging and determine, therefore, the state of 
each atom that leaves the cavity. 

In the successful case, furthermore, the next atomic 
pair is transported into the cavity and the next purifi- 
cation round takes place with the same stationary atom 
(permanent qubit). In the unsuccessful case, however, 
the stationary atoms have to be reinitialized and the en- 
tire sequence from Fig. [ija) should be restarted. The 
unsuccessful purification event, therefore, might demand 



a high amount of atomic pairs in the chain and require 
unreasonable operational times of the entire purification 
protocol. In order to avoid this potential bottleneck and 
perform the entanglement purification in a resource- and 
time-efficient way, in our previous, paper we suggested a 
procedure that exploits both (i) the periodicity of cavity- 
mediated evolution and (ii) non-destructive state mea- 
surements of atoms in the chain (see Sec. IV). 



III. HIGH-FIDELITY DYNAMICAL 
ENTANGLEMENT PURIFICATION 

As seen from Fig. |2|a), the output fidelity ^3(7,^2) 
(dashed curve) obtained after three successful purifica- 
tion rounds is still far from unit fidelity as required by a 
realistic quantum repeater. In fact, this output fidelity 
enables one to perform only a few swapping operations 
between the purified entangled pairs of neighboring re- 
peater segments until the fidelity of the resulting pair 
(distributed over a larger distance) drops to the initial 
fidelity /. Another bottleneck in our scheme is the neces- 



6 



sity to transfer the electronic population from the qubit- 
storage states to the cavity-active states (and backwards) 
in order to control the cavity- mediated evolution of atoms 
inside the cavity and protect our qubits against the de- 
coherence caused by the fast-decaying excited state |e) 
[see Fig. [TJc)]. Obviously, these two obstacles make our 
scheme less attractive to be considered in practice. 

In this section, we propose a modified purification 
scheme, in which we significantly improve the output fi- 
delity of remotely entangled atoms and get rid of the su- 
perfiuous laser pulses required to transfer the electronic 
population of atoms. By introducing one additional en- 
tanglement protocol in each repeater node and by op- 
timizing the laser beams required to control the entire 
scheme, we achieve an almost unit output fidelity after 



the same number of successful purification rounds. This 
dramatic improvement, therefore, allows for multiple en- 
tanglement swapping operations on the purified pairs. 

Similar to the original scheme that we presented in 
the previous section, the modified scheme includes two 
repeater nodes A and B sharing one permanent qubit 
pair Pab^ characterized by the density operator pj/, and 
a finite set of temporary entangled pairs as displayed in 
Fig. [sj^a). Each temporary entangled pair is given by the 
rank- two mixed state ([T]) in the basis {|0), |1)}, such that 
the fidelity ([2| of each pair is above the threshold value 
of 1/2. In contrast to the original scheme, however, right 
before the permanent pair is supplemented by the tempo- 
rary pairs Iab and 2ab^ these two (separate) entangled 
pairs are merged into the four-qubit entangled state 



Pf — 2 {\^1A,2A^^1B,2b){^1A,2A^^1B,2b\ + I'^1A,2A' '^1S,2b) ('^1A,2A' '^1S,2B l) 

^ 2(l-2/ + 2/2) (|^1A,2A'^1B,2b)(V^1A,2A'V^1B,2bI + \^1A,2A^ ^IB ,2b) {'I^1A,2A^ "I^IB ,2B l) ■ (12) 



This entangled state is generated using an additional 
entanglement protocol that occurs prior to the purifi- 
cation gate in our scheme. By this protocol, the pairs 
1^1,2^ and 1b, 2^ interact locally within the repeater 



nodes A and respectively, such that the state (12) 
is generated. In Fig. [sj^d) we display the experimental 
setup of our modified scheme. In contrast to the setup 
displayed in Fig. [TJd), we added (i) high-finesse cavities 
Ci and C3, (ii) photon detectors Di and and (iii) 
laser beams L2 and L5 to each of the repeater nodes 
A and B, respectively. These ingredients are compati- 
ble with the resources utilized in the original purification 
scheme and they form together the fusion block that is 
framed by a rectangle in Fig. [sj^d). Finally, the laser 
beams Li and L4 act continuously along the z-axis and 
together with the cavity field of Ci (C2) and C3 (C4), re- 
spectively, produce the two-photon (Raman) transition 
between the states |0) and |1) of the coupled atoms, such 
that the fast-decaying excited state |e) remains almost 
unpopulated (see below). 

Being transported from the entanglement distribution 
block into the cavity Ci (C3), the atoms 1^,2^ (1b, 2^) 
couple simultaneously to the same cavity mode and both 
laser beams L2 {L^) and Li (L4) as displayed in Fig.jsj^b). 
Assuming the non-zero cavity relaxation rate associ- 
ated with Ci (C3), the evolution of the coupled atom- 
cavity-laser system is governed by the master equation 
|29j 

p = — ^ [Hs^ p] -\- (2 a pa^ — ao) p — po)a) = C p] (13) 

Hs = ^{a + a^)[a^+a^), (14) 
where p is the density operator describing the state of the 



two atoms together with the cavity mode, C is the Lind- 
bladian superoperator that acts on the density operator, 
erf is the respective Pauli operator in the basis {|0), |1)}, 
and J2 is the coupling between the atoms inserted into 
the same cavity mode and subjected to the two laser 
beams. We show in Appendix [A] that the above Hamilto- 
nian is produced deterministically in our setup assuming 
both (i) the strong driving regime of atoms and (ii) the 
large detuning limit for laser and cavity fields. 

The evolution of atoms 1a, 2^ (1b, 2b) coupled to the 
field of cavity Ci (C3) due to Eq. ([l3| is completely de- 
termined by the exponent e^^ and the initial state of the 
cavity and the atoms. Since the pairs of atoms I^b and 
2ab are initially entangled, we have to consider the com- 
posite density operator 



PAB 



(p}0P^0|Oa,Ob)(Oa,Ob|) , (15) 



describing the state of two atomic pairs and two initially 
empty cavities at a given time t. In this expression, pj. 
and are the density operators of the entangled pairs 
l^B and 2ab, respectively, while |0a) and |0b) denote 
the vacuum states of cavities Ci and (^3, respectively. In 
Appendix [B] we show, moreover, that conditioned upon 
the no-photon measurement of the leaked cavity field in 
both repeater nodes, the state (15) reduces to the state 
(12) in the steady-state regime ^ 1), that is 



1,2 

pj 



(6^,6bIp^^b|Oa,Ob) 
1^[{QaM~Pab\^a,0b)] ' 



(16) 



where p^^^ is the operator ( 15 ) in the steady-state regime. 
The measurement of the leaked cavity field is performed 
using the photon detector Di {D2) that is connected to 
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FIG. 4. (Color online) (a) Fidelities /) (solid curve) and 
Fz^f^F'i) (dashed curve) given by Eqs. (20) and (21), respec- 
tively. The dotted curve displays the fidelity given by Eq. ([sf 
obtained in the original scheme, (b) Plot of F{f,n) in the 
modified scheme as a function of input fidelity / and number 
n of purification rounds, (c) Off-diagonal contributions Go{f) 
(dashed curve) and Gsif, G2) (solid curve) given by Eqs. (30) 
and (34), respectively. 



the neighboring repeater node through a classical com- 
munication channel. We stress that since the detection 
of the leaked cavity field in our scheme discriminates be- 
tween a vacuum state (no clicks) and a strong coherent 
state (many clicks), the efficiency of the detectors Di and 
D2 can take rather moderate values. 



Assuming that the four-qubit entangled state (12) has 
been successfully generated, the permanent pair Pab is 
supplemented by two temporary pairs Iab and 2ab as 
displayed in Fig. [sja). Similar to the original scheme, 
each repeater node contains now one triplet of qubits 
and each of these triplets evolves due to the isotropic 



Heisenberg XY Hamiltonian 
HxY = ^2^[^i 



cr, a, 



+1; ' 



i=l 



over the time period (n = 0, 1, 2, . . .) 



T=-,n- 



7-1 
^3 ' 



(17) 



(18) 



such that = cr^ and cf^ = , and where J3 is the 
coupling between the qubits. In Appendix [C| we show 
that the Hamiltonian ( pT] ) is produced deterministically 
in our scheme by coupling simultaneously three atoms to 
the same cavity mode C2 {C4) and the laser beam Li 
(L4) in the large detuning limit [see Fig. [sjc)]. 

Similar to the original scheme, this evolution is fol- 
lowed by the projective measurement of qubit pairs Iab 
and 2ab in the basis {|0), |1)} and the exchange of the 
projection outcomes between the two repeater nodes by 
means of classical communication. The entanglement pu- 
rification is successful if the outcome of projections agrees 
with ([5|. In this case, the (unprojected) permanent qubit 
pair is described again by the density operator ([6|, where 

F^fj')- [25-50/ + IMP) f 



169 + 194/2 - 144/' + (288/' - 338)/ ' 

(19) 

such that Fipp) = F{f,f) > F{pf,). In Fig.ga), we 

compare the fidelity 



F{fJ) 



(25 -50/ + 194/2)/ 
169 - 482/ + 482/2 



(20) 



(solid curve) with the respective fidelity given by Eq. ([8| 
(dotted curve). We see that the growth of fidelity in the 
modified scheme is almost twice as large as in the origi- 
nal purification scheme. This nice result, however, relies 



merely on the input state (12) that is entangled strongly 
if compared to the separable state pj- (g) p'j used in the 
original scheme. We recall that our scheme relies on the 
effect of entanglement transfer between the networks of 
evolving spin chains introduced in Ref. [27] and realized 
in our scheme using the cavity QED framework. The 
stronger the entangled state that we provide as the in- 
put for the purification block in our scheme, the more 
entanglement is transferred to the permanent qubit pair. 

Assuming that each purification round is successful, 
the sequence from Fig. [sj^a) leads to a gradual growth of 
entanglement fidelity (of stationary atoms) with regard 
to the respective fidelity obtained in the previous round. 
Similar to the original scheme, we analyze quantitatively 
the sequence (10) in order to understand how much the 



output fidelity increases with each purification round. In 
Fig. gb), we show the plot of F(/, n) that describes 
the difference between the final fidelity F^(/, F^_i) ob- 
tained after n (successful) rounds and the initial fidelity / 
(n = 0). For / > 0.75, this function exhibits a dramatic 
growth during the first three rounds that saturates and, 
with increasing n, yields a negligible growth with regard 
to the following fixed point fidelity 



f (25 - 50/ + 194/2 



4826809 - 33772038/ + 103411314/2 - 179097440/^ + 189095940/* - 119456664/5 + 39818888/^ 



(21) 



displayed in Fig. [4]^a) by a dashed curve. We readily see 
that we obtain an almost unit output fidelity for / > 0.75 
after the same numner of purification rounds as in the 
original scheme. For / < 0.75, however, the function 
n) continues to grow with each purification round. 
In this case, the optimal number of purification rounds 
and the respective fixed point fidelity has to be deter- 
mined for each particular value of / separately. 



A. Evolution governed by Hamiltonian (17) and 
the purification gate 

In this section, we already explained that the density 
operator (|6| with the function (19) characterize com- 



pletely the permanent qubit pair obtained in our (mod- 
ified) scheme after a single purification round. In this 
subsection, we analyze briefiy the evolution governed by 
the Hamiltonian (17) and connect it with the main re- 



sults utilized in this paper. 



The atomic evolution governed by Hamiltonian (17) 



/c=l 



\k){k\, 



(22) 



is completely determined by the energies E}^ and vec- 
tors which satisfy the eigenvalue equality Hxy \k) = 
Ej^ \ k) with orthogonality and completeness relations 
{k\k^) = Skk' and ^\k){k\ = /, respectively. With the 
help of Jordan- Wigner transformation [30], this eigen- 
value problem can be solved exactly (see, for instance, 
Ref. [31]). Since the evolution operator ( [22] ) acts on the 
states of one atomic triplet that is entangled with an- 
other atomic triplet in the neighboring node, we have to 
consider the composite evolution operator 



U{t) 



E 

k,k'=l 



Earlier we explained that right before each atomic pair 
from node A enters the cavity, it becomes entangled with 
another atomic pair from node B, such that the four- 
qubit state (12) is generated. We denote the density 



operator of the stationary atoms by pj,. According 



to 



the evolution operator (23) and this notation, the state 



of both atomic triplets in nodes A and B is described by 
the six-qubit density operator 



^1,2,3 



{t,f,f) = U{t)[py^p},)u\t) (24) 



nodes A and B is described by the density operator 
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pi'^'^(T,/,/) = ^ p]f^'{TJJ') |v,)(v,|, (25) 



where 2^ composite vectors |v^), satisfying the orthog- 
onality and completeness relations (v^lvj) = 5ij and 
— ^5 respectively, have been introduced. 
In order to finalize one purification round, the conveyed 
atomic pairs are projectively measured, such that the 
projected density operator 



p(T, /,/')= 



a, 13=1 



p'J'HtJJ') 

-PsUCc(^7 fj /O 



|va)(v^|, (26) 



with 



Psucc(T,/,/') = Tr 



a, 13=1 



describes the state of the stationary atoms. In the above 
expressions, the Greek indices run over the four different 
values given by 

|v«) = (0iA,l2A,0iB,l2B|v«) 7^0, (27a) 
or |v«) = (1ia,02a,1ib,02b|v«) 7^0, (27b) 

which correspond to the outcomes of the projections ([5|. 

Using the six-qubit density operator (24), we have 
routinely computed the matrix elements p-j '^(T^f^f) 
which, however, are rather bulky to be displayed here. 
With the help of these matrix elements, we confirmed 
that the density operator (26) coincides with the rank- 
two mixed state (|6), where the function F(/, /') is given 
by Eq. (19). Another relevant result that we learn from 



the expression (|26|) is the cyclic property (m = 1, 2, . . .) 



p{m7rJ^\fJ') = p{0JJ')=p}, 



(28) 



In the case of an unsuccessful purification event, we ex- 
plain in Sec. |IV| how this property makes it possible to re- 
store the output fidelity to the value obtained in the pre- 
vious (presumably successful) purification round. This 
restoration enables one to counteract unsuccessful purifi- 
cation events and to implement our scheme in a resource- 
and time-efficient way. 



B. Remarks on the entanglement distribution 
between stationary atomic qubits 



that evolves over the time period T given by Eq. (18). 
After this evolution, the state of both atomic triplets in 



Throughout the paper, we assumed that the station- 
ary atoms are initially entangled, such that the fidelity is 
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above the threshold value of 1/2. In our previous paper, 
we suggested that there is no need to introduce an ad- 
ditional entanglement distribution protocol in our setup 
in order to entangle the stationary atoms prior to the 
purification. Instead, it was suggested to prepare ini- 
tially both permanent atoms in the ground state and run 
our purification scheme. We showed that one success- 
ful purification round induces the entanglement of the 
stationary atoms and ensures that the fidelity of result- 
ing density operator (almost) coincides with the fidelity 
/ of the temporary pairs and pj. In other words, 
one purification round entangles two (initially separable) 
stationary atoms, such that the fidelity of the temporary 
pairs is mapped to the fidelity of the stationary density 
operator. 

In our modified scheme, we utilize the same procedure 
as in the original scheme. We generate two (separate) 
entangled pairs and let them be conveyed through the 
cavity Ci (C3), such that the (probabilistic) entangle- 
ment protocol that produces the four-qubit state (12) is 



which are motivated by the optimal number of purifica- 



switched off. It can be shown that a successful purifica- 
tion round with these two entangled pairs transforms the 
state of the permanent atoms (prepared initially in the 
ground state) into an entangled state described by 



G, 



0{f){\<l>XB){<pA,B\ + \<pA,B){<l>XB\) . (29) 



with 



Foif) = 



1 + 48/ + 32/2 



Goif) = 



9 - 32/ + 32/2 



(30) 

where cq = 82, ci = C2 = 64 are the only non-zero coeffi- 
cients. 

Obviously, the above state is no longer a rank- two 
mixed state like in Eq. ([T]) because of the off-diagonal con- 
tribution Gq(J) displayed in Fig.CTc) by a dashed curve. 
The fidelity Q associated with ( [29| ), however, is slightly 
larger than the fidelity / associated with the temporary 
pairs and pj. The role of this initialization round is 
solely to entangle the stationary atoms and, therefore, 
it has to be followed by a number of purification rounds 
leading to the gradual growth of fidelity, 

FoU) < Fi{f,Fo) < ... < F„(/,F„_i), (31) 

and the gradual reduction of (off-diagonal) contributions, 

Goif) > Gi(/,Go) > . . . > Gn(/,Gn-i) . (32) 

Using the above sequences, we calculated the fidelity 

(1 + 48/ + 32/2)(25 - 50/ + 194/2)3 

and the respective off-diagonal contribution 

274625 (9 - 32/ + 32/2) (1 _ 2/ + 2/2)3 



(33) 



G3(/,G2) = 



tion rounds obtained previously [see (11) and (21)], and 
where 

^^0 = 195493577; di = 1442887766; ^2 = 4716352898; 
ds = 8883640864; d^ = 10517241220; d^ = 7944708952; 
de = 3738576328; dr = 934577152; ds = 233644288 

are the only non-zero coefficients. 

In contrast to the dashed curve describing Go(/) 
in Fig. [4]^c), the solid curve describing (34) deviates 
slightly around the constant value of 0.004. To a good 
approximation, therefore, the off-diagonal contribution 
G3(/, G2) can be neglected and the resulting density op- 
erator takes the form of the rank-two mixed state ([T|). 
We have verified, moreover, that the output fidelity ([33]) 
(almost) coincides with the output fidelity (21) displayed 
in Fig.pja) by a dashed curve. The price we pay for one 
extra (successful) purification round prior to the main 
sequence of rounds, therefore, is clearly compensated by 
the more moderate demand of physical resources in our 
purification scheme. 



IV. COUNTERACTING THE UNSUCCESSFUL 
PURIFICATION EVENTS 

Throughout the paper, we assumed that each purifica- 
tion round is finalized successfully leading to a gradual 
growth ([9| or ( [31] ) of entanglement fidelity. In the case 
of an unsuccessful purification event, i.e., when the out- 
comes of the projective measurement disagree with ([5|, 
the stationary atoms should be reinitialized and the en- 
tire scheme restarted. Since the probability to get the 
above two (out of 2^) combinations of projective measure- 
ments is rather small, multiple unsuccessful events might 
require a large amount of atomic pairs in the chain and 
unreasonable operational times. In order to avoid this 
bottleneck and perform the entanglement purification in 
a resource- and time-efficient way, in our previous paper 
we suggested a procedure that exploits (i) the periodicity 
of the cavity-mediated evolution and (ii) non-destructive 
state measurements of the atoms in the conveyed chain. 

Our procedure is as follows. We set the velocity of the 
atomic chain that provides the temporary qubit pairs, 
such that the transition time of these atoms through 
the cavity is higher than the operational time (18) re- 



(34) 



quired for one single purification gate. In this case, an 
atomic pair is still located inside the cavity after it under- 
goes (together with the stationary atom) one purification 
gate. Now instead of waiting until this pair leaves the 
cavity and is project ively measured outside, we probe 
the atomic states inside the cavity by means of non- 
destructive state measurements, which are based on the 
cavity transmission technique (see below). By exchang- 
ing the outcomes of these measurements between both 
repeater nodes, the state of stationary and temporary 
atoms (which are still in the cavity) can be restored to its 
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initial product state in the case of an unsuccessful purifi- 
cation event using the periodicity of the cavity-mediated 
evolution (28). 



More specifically, the restoration is realized by switch- 
ing on the cavity- mediated evolution (17) once more after 



the purification gate and the non-destructive measure- 
ments. According to Eq. (28), the operational time of 



this additional evolution is set to {mirJ^^—T), such that 
the quantum state of the stationary atomic pair evolves 
to its initial state p^,. This auxiliary evolution restores 
the initial state of the stationary atoms in each repeater 
node and, as a consequence, the entanglement fidelity 
that has been obtained in the previous (presumably suc- 
cessful) purification round. With the help of such non- 
destructive measurements and one extra cavity-mediated 
evolution, therefore, the entanglement fidelity of the sta- 
tionary atoms either increases due to successful purifica- 
tion or remains unchanged due to unsuccessful purifica- 
tion. 

In our experimental setup, the non-destructive mea- 
surements of the atomic states in each repeater node are 
performed with the help of (i) a probe laser beam that 
acts along the y-ax.is of cavity C2 {C4) and is tuned to the 
cavity resonance frequency and (ii) a detector that mon- 
itors the transmission of this probe beam. The physical 
mechanism behind these non-destructive measurements 
exploits the suppression of cavity transmission that arises 
due to the strong atom-cavity coupling (see, for instance, 
Refs. [SH [33]). Recall that each atom in our scheme 
is a three-level atom in the A-type configuration [see 
Figs. [3]^b),(c)], where only the states |0) and |e) are cou- 
pled to the cavity field. Apparently, if one such atom 
prepared in the state |0) is inserted into a cavity that 
is sufficiently detuned from the atomic |0) ^ |e) transi- 
tion frequency, then the cavity transmission drops in de- 
pendence on the atom-cavity detuning and atom-cavity 
coupling. On the other hand, the cavity transmission re- 
mains unaffected if an atom prepared in the state |1) is 
inserted into the cavity. 

Once sufficiently many readouts of the cavity trans- 
mission are recorded, the mechanism described above en- 
ables one to determine (non-destructively) the state of a 
single atom with a good efficiency [32]. Since the atom- 
cavity coupling increases proportionally with the number 
of inserted atoms, the same mechanism enables one also 
to distinguish the following three composite states: (i) 
|0i,02), (ii) |0i,l2) or |li,02), and (iii) |li,l2) of the 
two (temporary) atoms which are transported and cou- 
pled to the cavity C2 {C4) [33 . An alternative (non- 
destructive) detection scheme that fits nicely to our ex- 
perimental setup was proposed in Ref. [34]. Similar to 
the previous approach, this scheme utilizes three-level 
atom(s) which are inserted into a detuned optical res- 
onator, such that the atomic state |1) is decoupled from 
the cavity field. In this scheme, however, the field pro- 
duced by the probe laser beam is transmitted through 
the coupled atom-cavity system and measured by means 
of a homodyne detection. This homodyne detection re- 



veals the changes of phase of the transmitted probe beam, 
where the phase shift is proportional to the number of 
atoms in the state |0). This scheme, therefore, enables 
one to distinguish the above three composite states (i), 
(ii), and (iii) of the two atoms which are coupled to the 
cavity. 

The approach presented in this section requires that 
short atomic chains are transported with a constant ve- 
locity along the experimental setup and coupled to the 
cavity-laser fields in a well controllable fashion. For this 
purpose, we introduced in our setup [see Fig. [sj^d)] (i) 
a magneto-optical trap (MOT) that plays the role of an 
atomic source and (ii) an optical lattice (conveyor belt) 
that transports atoms into the cavity from the MOT with 
a position and velocity control over the atomic motion. 
The proposed setup is compatible with existing exper- 
imental setups [351137] . in which the above devices (i) 
and (ii) are integrated into the same framework together 
with a high-finesse optical cavity. The number-locked in- 
sertion technique [38], moreover, enables one to extract 
atoms from the MOT and insert a predefined pattern of 
them into an optical lattice with a single-site precision. 
It was already demonstrated that an optical lattice pre- 
serves the coherence of transported atoms and can be 
utilized as a holder of a quantum register. By encoding 
the qubits by means of hyperfine atomic levels, a qubit 
storage time of the order of seconds has been demon- 
strated within this register [28l [39] . 



V. SUMMARY AND DISCUSSION 

In this paper, an efficient, high-fidelity scheme was pro- 
posed to purify the low-fidelity entangled atoms trapped 
in two remote optical cavities. This scheme is a modi- 
fication of the purification scheme proposed in our pre- 
vious paper [24] that exploits the natural evolution of 
spin chains instead of CNOT gates. Similar to the orig- 
inal scheme, the modified scheme uses a cavity- QED 
framework, namely (i) short chains of atoms, (ii) high- 
finesse optical cavities, and (iii) detectors for the projec- 
tive measurement of atomic states. In contrast to the 
original scheme, however, one additional entanglement 
protocol was introduced in each repeater node, and the 
laser beams which are used to control the entire scheme 
were optimized. With the help of these modifications, an 
almost unit output fidelity was achieved after the same 
number of successful purification rounds as in the original 
scheme. 

Similar to the original paper, furthermore, the mod- 
ified scheme was supplied with a detailed experimental 
setup, and a complete description of all necessary steps 
and manipulations has been given. A comprehensive 
analysis of fidelities obtained after multiple purification 
rounds was performed and the optimal number of rounds 
was determined. An additional mechanism was suggested 
that enables one to counteract unsuccessful purification 
events and ensures that our scheme is implemented in 
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a resource- and time-efficient way. Following recent de- 
velopments in cavity QED, moreover, we briefly pointed 
to and discussed through the paper a few practical is- 
sues related to the implementation of our purification 
scheme. We stress that, although the proposed purifi- 
cation scheme is experimentally feasible, its complete 
experimental realization (including the non-destructive 
measurements) still poses a big challenge. 

The high-fidelity purification scheme proposed in this 
paper enables one to perform multiple entanglement 
swapping operations and thus opens a route towards an 
efficient and experimentally feasible quantum repeater 
for long-distance quantum communication. More specif- 
ically, in our experimental setup, each atom in node A 
has to be entangled with another atom from node B right 
before they enter the cavities Ci and C3 for further pro- 
cessing. The (low- fidelity) entanglement between these 
atoms is distributed non-locally using the entanglement 
distribution block indicated in Fig. IsFd) by a rectangle. 
In order to entangle two (three-level) atoms located at 
distant repeater nodes A and B, we find the entangle- 
ment distribution scheme proposed in Ref. [40 the most 
appropriate. This scheme is also realizable in the frame- 
work of cavity-QED and, therefore, it utilizes the same 
physical resources as our purification scheme. 

By this scheme, a coherent-state light pulse interacts 
with the coupled atom-cavity system in node A, such 
that the optical field accumulates a phase conditioned 
upon the atomic state in this node. Afterwards, the light 
pulse propagates to node B, where it interacts with the 
second coupled atom-cavity system and accumulates an- 
other phase conditioned upon the atomic state in this 
node. The resulting density operator [25 



with 



A,B/ 



1 



\ + {l-mA.B){4>A.B\ 



(35) 



± ^ e-'^^Ci)\lA,OB) + \ e"'1C2>|0^, Is) , 

describes the state of both atoms and the coherent light 
pulse, where Co = ^J^OL^ C\ = ^ote^^ ^ C2 = ^Jfjae~''^ 
denote the phase-rotated and channel-damped coherent 
state a, e = sin (9, while / = (1 + e-(i-^)^'(i-^°«^))/2 
plays the role of the entanglement fidelity. 

The resulting (phase-rotated) coherent pulse becomes 
disentangled from the atoms with the help of homo- 
dyne detection followed by post-selection [^ or, alterna- 
tively, using unambiguous state discrimination [2 53 . This 
projects the state (35) onto an entangled state of two 
atoms that coincides with the rank- two mixed state ([T]). 
We remark that the conditioned phase rotation exploited 
in the entanglement distribution scheme is naturally re- 
alized in a cavity-QED framework using the single atom- 
cavity evolution in the dispersive interaction regime. 
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Appendix A: Derivation of the Hamiltonian (14) 



In this appendix, we show that the Hamiltonian (14) 



is produced deterministically in our setup. Specifically, 
two (three-level) atoms are subjected to the field of the 
(initially empty) cavity Ci (C3) and the fields of laser 
beams Li (L4) and L2 {L^) simultaneously as displayed 
in Fig. [sj^b). The evolution of this coupled atom-cavity- 
laser system is governed by the Hamiltonian (/c = 1, 2) 



Hi = hujca"^ a - [|«|e)fe(0| 



(Al) 



-iojp t 



e)k{0\)-H.c, 



+ [cji|l)fc(l| +a;£;|e)fc(e| +cJo|0)fc(0|] , 

k 

where g denotes the coupling strength of an atom to the 
cavity mode, while ft denotes the coupling strengths of 
an atom to both laser fields. 

We assume that ujc = ojp and switch to the interaction 
picture using the unitary transformation 

IJ_^ ^ g-a[E(c^i|l>fc(l|+c^s|e)fc(e| + (u;z.+u;i-u;p)|0)fc(0|)+u;pa"t"a] 

We assume, moreover, that Al = Ac = —A, where the 
notation Al = {uje —^1) — and Ac = {(jOe —^0) —^c 
has been introduced. In t he a bove interaction picture, 
therefore, the Hamiltonian (Al) takes the following form 



Ho 



2' 



k 

-iAt 



2 



a\e)k{ 



(A2) 



{\e)k{l\^\e)km-H.c. 



We require that A is sufficiently far detuned, such that 
no atomic |e) ^ |0) and |e) ^ |1) transitions can occur. 
We expand the evolution governed by the Hamiltonian 
(A2) in series and keep the terms up to the second order. 



Uo 



I 



H2dt' 




H2 H2dt" \ dt' . 



By performing integration and retaining only linear-in- 
time contributions, we express this evolution in the form 



U2 



I 



exp 



(A3) 



where the effective Hamiltonian is given by 



g\l)kma + H.c.] , (A4) 
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after removing the constant contributions. 

At this stage, we switch from the atomic basis {|0), |1)} 
to the basis {|+), | — )}, where 



W = ^(|0) 



-> = ^(|o> 



(A5) 



In this basis, the Hamiltonian (|A4|) takes the form 



(A6) 



+g (sUa + a^) + {Sl-Sk){a-a))) 

where Sk = \-)k{+\ and 5"^ = |+)fc(+| - |-)fe(-|, and 
where we removed any constant contributions. We switch 
once more to the interaction picture with respect to the 
first term of (|A6|. In this interaction picture, we obtain 



at) 



8A 



k 



5fee-^fe*)(a-at) 



(A7) 



In the strong driving regime, i.e., for Vt ^ {^, A}, we 
ehminate the last (fast oscihating) term using the same 
arguments as for the rotating wave approximation. The 
Hamiltonian (A7), therefore, reduces to 



(A8) 



where we used the identity = cj^ . The resulting 



Hamiltonian (A8) coincides with the Hamiltonian (14) 



under the notation J2 = gVt/{A/S). 



Appendix B: Steady-state solution of Eq. (13) 



In this appendix, we show that the mixed state ( 12| is 
conditionally generated by means of evolution (13) that 



takes place simultaneously in both repeater nodesfA and 
B) with the initial state 



|0a,0b)(0a,0b| 



(Bl) 



in the steady-state regime. In Sec. HI, we considered the 
expression ([l5| based on the exponent e^^ that is difficult 
to evaluate. Here we apply sequentially the steady-state 
solution of Eq. (13) to the coupled (atom-atom-cavity) 
system in node A and, afterwards, to the coupled system 
in node B. 

In order to proceed, we consider ffist the solution of 



(13) for an initially empty cavity Ci and a general mixed 
state of two atoms in node A (i, j = 1, . . . , 4) 



(p°),,|uf){u/|®|0^)(0^| 



(B2) 



where we switched to the atomic basis (A5), such 
that \uf) = |+ia,+2a), |u^) = \-iA,-2Ah\ui) = 



l+iA, -2a), and |uf ) = \-ia, +2a) form together an or- 
thogonal basis. To our best knowledge, the master equa- 
tion (13) was solved in Refs. [4TJ [42] only for an initial 
pure state of atoms, which is not appropriate to be used 
in our case. We therefore (re) solved this master equation 
for an initial mixed state of two atoms (B2). Assuming 



the strong atom-cavity coupling J2 ^ the solution we 
found in the steady-state regime Kt ^ 1 can be expressed 
as 

Pa = ^{PAh A., \uf){uf\0\-ui a^J){-u, <| , (B3) 
where = 2LJ2/tv is the amplitude of the coherent 



state, 
U4 = 0. 



^^.)^o with ui = 1, U2 = -1, and us 



Recall that the evolution (13) takes place simultane- 



ously in both repeater nodes A and B with the initial 
state (|B1|), which we cast into the form 

f^AB 



E (pABhki |uf ,uf)(u/,uf I ® |0^,0b)(0^,0b| 



i,j,k,l 

Y,{p'^),,\uf){uf\^\OA){OA\ 



where we introduced the matrices 

(pAshkl = (uf , uf I ® p} \uf, uf ) 

)ijkl 



(B4) 



(B5) 



iPAh = J2(PABhki |uf ){uf I ® IObXObI . (B6) 



By identifying the expressions (B2) and (B4), we con- 



clude that the density operator (B3) gives the steady- 



state solution of (13) obtained in node A for the initial 



state ( |B3[ ). Conditioned upon the no-photon measure- 
ment of the leaked field from the cavity Ci, this steady- 
state solution reduces to 

{OaIp^IIOa) = J2(pAh K e-l"^"!'-"^ |u,^)(u/| 

^E(p^)fcHuf)(uf|®|OB)(OB|, (B7) 



where 

{p%)ki 



y^XPAB)ijkl \j e 



uf)(u 



(B8) 



Similar to node A, the density operator 

Pb = Y.^P%)ij K Iuf )(uf I c^B){-Uj <l . 



(B9) 

gives the steady-state solution of (13) for node B with 
the initial state (B8). Conditioned upon the no-photon 



measurement of the leaked field from the cavity C3, this 
solution reduces to 

(0a,0b|p^^b|0a,0b) 

k,l 

= E (p°A).,.K-l"^^l''-'''|u.^uf)(u/,uf |,(B10) 

i,j,k,l 



13 



where the notation Oij^i = UiUj + uj^ui^ {p^^)ijki = 
{PAB)ijki ^ij^kh and = = a^^ was introduced. 
The atom-cavity strong couphng ensures that ^ 1. 

To a good approximation, this observation impUes that 



the exponent in (BIO) vanishes for Oijki ^ 



(0a,0b|p^^b|0a,0b) 



^£ (pOA),,,,|uf,uf)(u/,uf|. 

(Bll) 

Owing to the exphcit form of (Bl), we have routinely 
computed the matrix elements {pABJijki which, however, 
are rather bulky to be displayed here. With the help 
of these matrix elements, we obtained the normalized 
density operator (Bll) in the form 



1,2 



2 - 4/ + 4/2 



(|'^^4,1B''^^4,2b)('^^4 ,1B'^2A,2bI + \^1A,1B^^2A,2b){^1A,1B^^2/ 



2A,2B\ 



'I^IAAB^ '^2A,2b)(^1A, 



'2A,2bI 



' ^2A,2b)('^LA,1B' '^2A,2b|) 
2 _ 4j ^ 4J2 ^2A,2b)(^1A,1B' ^2A,2bI + \'^1A,1B ^ '^2A,2b) i'^lA^lB ^ ^2A,2b\ 

~\^lAAB-> ^2A,2b){^1AAB^ ^2A,2B AB-> ^2A,2b){^IAAB-> ^2A,2B 



(B12) 



that describes a four-qubit entangled state and coincides 
with the state (12). 



Appendix C: Derivation of the Hamiltonian (17) 



In this appendix, we show that the Hamiltonian (17) 



is produced deterministically in our setup. Specifically, 
three (three-level) atoms are subject to the field of the 
(initially empty) cavity C2 (C4) and the field of laser 
beam Li (L4) simultaneously as displayed in Fig. [sj^c). 
The evolution of this coupled atom-cavity-laser system is 
governed by the Hamiltonian (k = 1, 2, 3) 



Hi — hujc a 



-iooL t 



e)k{l\-H.c. 



LhY^ |a|e)fe(_, , 2 

k L 

^ujE\e)k{e\ +a;o|0)fe(0|] , 



where the notation Al = {oje — ^1) — and Ac = 
{ue — ^0) — has been introduced. 

We require that and Ac are sufficiently far de- 
tuned, such that no atomic |e) ^ |0) and |e) |1) tran- 
sitions can occur. W e expand the evolution governed by 
the Hamiltonian ( C2 ) in series up to the second order. By 



performing integration and retaining only linear-in-time 
contributions, we express this evolution in the form (A3), 



where the effective Hamiltonian is given by (we assume 
that the cavity field is initially in the vacuum state) 



Hs = hAa^ a- 



H.c], (C3) 



(CI) where A = A^ — Ac- We switch one more time to the in- 



teraction picture with respect to the first term of ( C3 ). In 



this interaction picture, the Hamiltonian takes the form 



4Ar ^ L 



iAt\ 



l)k{0\^ H.c]. (C4) 



where g denotes the coupling strength of an atom to the 
cavity mode, while Q denotes the coupling strength of an 
atom to the laser field. 

We switch to the interaction picture using the unitary 
transformation 

^ g-a[E(c^i|l>fc(l|+c^s|e)fc(e|+u;o|0)fc(0|) + (u;i+u;i.-u;o)a'^a] 



In this picture, the Hamiltonian (CI) takes the form 
H2 = h (Al - Ac) a'^a 



(C2) 



jAr t 



e)k{0\ 



iArt 



e)k{l\-H.c. 



We require, finally, that A is sufficiently far detuned. 
As above, we expand again the evolution governed by 
the Hamiltonian (C4) in series up to the second or- 



der and retain only linear-in-time contributions after 
the integration. This leads to the effective Hamiltonian 
(^,j = 1,2,3) 



16A| A 



H. 



^|Oi,li)(li,0,-| + ^|l)fc{l| 



(C5) 



Since the second term in this Hamiltonian commutes 
with the ffrst term, we eliminate the second term by 
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means of an appropriate interaction picture. The re- 
sulting Hamiltonian, i.e., the first term of (C5), co- 



incides with the Hamiltonian (17) under the notation 



[1] W. K. Wootters and W. H. Zurek, Nature 299, 802 
(1982). 

[2] D. Dieks, Phys. Lett. A 92 A, 271 (1982). 

[3] H.-J. Briegel, W. Diir, J. I. Cirac, and P. Zoller, Phys. 

Rev. Lett. 81, 5932 (1998). 
[4] C. H. Bennett, G. Brassard, S. Popescu, B. Schumacher, 

J. A. SmoUn, and W. K. Wootters, Phys. Rev. Lett. 76, 

722 (1996). 

[5] D. Deutsch, A. Ekert, R. Jozsa, C. Macchiavello, 
S. Popescu, and A. Sanpera, Phys. Rev. Lett. 77, 2818 
(1996). 

[6] M. Zukowski, A. Zeilinger, M. A. Home, and A. K. Ekert, 

Phys. Rev. Lett. 71, 4287 (1993). 
[7] Z. Zhao, T. Yang, Y.-A. Chen, A.-N. Zhang, and J.- 

W. Pan, Phys. Rev. Lett. 90, 207901 (2003). 
[8] R. Reich et al., Nature 443, 838 (2006). 
[9] T. Yang et al., Phys. Rev. Lett. 96, 110501 (2006). 
[10] H. de Riedmatten et al., Phys. Rev. A 71, 050302(R) 

(2005). 

[11] Z.-S. Yuan, Y.-A. Chen, B. Zhao, S. Chen, J. Schmied- 

mayer, and J.-W. Pan, Nature 454, 1098 (2008). 
[12] C.-W. Chou et al., Science 316, 1316 (2007). 
[13] Y.-B. Sheng, F.-G. Deng, and H.-Y. Zhou, Phys. Rev. A 

77, 042308 (2008). 
[14] B. Zhao, M. Miiller, K. Hammerer, and P. Zoller, Phys. 

Rev. A 81, 052329 (2010). 
[15] C. Wang, Y. Zhang, and G.-S. Jin, Phys. Rev. A 84, 

032307 (2011). 

[16] N. Sangouard, C. Simon, H. de Riedmatten, and 
N. Gisin, Rev. Mod. Phys. 83, 33 (2011). 

[17] P. van Loock, Laser Photonics Rev. 5, 167 (2011). 

[18] W. Diir, H.-J. Briegel, J. L Cirac, and P. Zoller, Phys. 
Rev. A 59, 169 (1999). 

[19] L. Isenhower et al., Phys. Rev. Lett. 104, 010503 (2010). 

[20] A. Gabris and G. S. Agarwal, Phys. Rev. A 71, 052316 
(2005). 

[21] S.-B. Zheng and G.-C. Guo, Phys. Rev. Lett. 85, 2392 
(2000). 



[22] T. Tanamoto, K. Maruyama, Y.-X. Liu, X. Hu, and 
F. Nori, Phys. Rev. A 78, 062313 (2008). 

[23] N. Schuch and J. Siewert, Phys. Rev. A 67, 032301 
(2003). 

[24] D. Gonta and P. van Loock, Phys. Rev. A 84, 042303 
(2011). 

[25] P. van Loock, N. Liitkenhaus, W. J. Munro, and 
K. Nemoto, Phys. Rev. A 78, 062319 (2008). 

[26] E. Lieb, T. Schultz, and D. Mattis, Ann. Phys. (N.Y.) 
16, 407 (1961). 

[27] A. Casaccino, S. Mancini, and S. Severini, Quant. Inf. 

Process. 10, 107 (2011). 
[28] S. Kuhr et al., Phys. Rev. Lett. 91, 213002 (2003). 
[29] P. Meystre and M. Sargent, Elements of Quantum Optics, 

(Springer Verlag Berlin Heidelberg, 2007). 
[30] P. Jordan and E. Wigner, Z. Phys. 47, 631 (1928). 
[31] X. Wang, Phys. Rev. A 64, 012313 (2001). 
[32] A. D. Boozer, A. Boca, R. Miller, T. E. Northup, and 

H. J. Kimble, Phys. Rev. Lett. 97, 083602 (2006). 
[33] M. Khudaverdyan, W. Alt, T. Kampschulte, S. Reick, 

A. Thobe, A. Widera, D. Meschede, Phys. Rev. Lett. 

103, 123006 (2009). 
[34] A. S. S0rensen and K. M0lmer, Phys. Rev. Lett. 91, 

097905 (2003). 

[35] S. Nussmann et al., Phys. Rev. Lett. 95, 173602 (2005). 
[36] K. M. Fortier, S. Y. Kim, M. J. Gibbons, P. Ahmadi, and 

M. S. Chapman, Phys. Rev. Lett. 98, 233601 (2007). 
[37] M. Khudaverdyan et al., New J. Phys. 10, 073023 (2008). 
[38] M. Karski et al., New J. of Phys. 12 065027 (2010). 
[39] D. Schrader, L Dotsenko, M. Khudaverdyan, Y. Mirosh- 

nychenko, A. Rauschenbeutel, and D. Meschede, Phys. 

Rev. Lett. 93, 150501 (2004). 
[40] P. van Loock et al., Phys. Rev. Lett. 96, 240501 (2006). 
[41] M. Bina, F. Casagrande, and A. LuUi, Opt. and Spec- 

trosc. 108, 356 (2010); Phys. Rev. A 77, 033839 (2008). 
[42] P. Lougovski, E. Solano, and H. Walther, Phys. Rev. A 

71, 013811 (2005); P. Lougovski, F. Casagrande, A. Lulh, 

and E. Solano, Phys. Rev. A 76, 033802 (2007). 



